Electron-phonon interaction and full counting statistics in molecular junctions 
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The full counting statistics of a molecular level weakly interacting with a local phonon mode is 
derived. We find an analytic formula that gives the behavior of arbitrary irreducible moments of 
the distribution upon phonon excitation. The underlying competition between quasi-elastic and 
inelastic processes results in the formation of domains in parameter space characterized by a given 
sign in the jump of the irreducible moments. In the limit of perfect transmission, the corresponding 
distribution is distorted from Gaussian statistics for electrons to Poissonian transfer of holes above 
the inelastic threshold. 
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It is now well established that the current noise gen- 
erated by electronic nanodevices contains valuable infor- 
mation on microscopic transport processes not available 
from measurements of the current- voltage characteristics 
A full characterization of the transport properties 
of a device requires not only the knowledge of current- 
current correlations, but rather the full counting statis- 
tics (FCS) has to be determined This amounts to 
determine the whole probability distribution P(„ (g) that 
a given charge q is transmitted through the device during 
a certain measurement time t^. 

Studies of FCS have been mainly restricted to non- 
interacting systems. Notable examples of such studies 
are single channel conductors Q or double quantum dot 
systems Q ■ The case of FCS in the presence of electron- 
electron interactions in the coherent transport regime has 
been much less analyzed In particular, the Kondo 
regime in quantum dots has been addressed in Ref.Q. 
On the other hand, molecular electronics is becoming 
a field of intense research activity In this case the 
coupling to vibrational modes plays an important role 
and provides an additional source for electronic correla- 
tions which may affect the counting statistics [1]. The 
case of atomic chains suspended between metallic elec- 
trodes provides another test system to analyze the ef- 
fects of electron-phonon coupling in transport properties 
0. At low temperatures (quantum regime) the onset of 
phonon emission processes is signaled by a brup t jumps 
in the system differential conductance [10|, llli- When 
certain conditions are met, the behavior of the conduc- 
tance jumps is entirely controlled by the transmission 
probability, evolving from a drop in conductance at high 
transmission to an increase at low transmission. These 
predictions were quantitatively confirmed in recent ex- 



periments [I4I . A natural question which arises concerns 



the behavior of noise and, more generally, the FCS for 
energies corresponding to the excitation of vibrational 
modes. Although some works have been devoted to the 
analysis of noise in the presence of e-ph coupling [l^, 
none of them tackled the problem of the determination 



of the FCS in the presence of e-ph interaction. 
The aim of this work is to study how the FCS of a 
molecular junction is modified by the coupling to a vibra- 
tional mode. On the basis of a simple model, we derive 
a compact analytical expression encoding the FCS for 
the experimentally relevant regime of weak e-ph interac- 
tions and strong coupling to the leads, which corresponds 
to the conditions of Ref. 12]. This expression allows to 
analyze the change of arbitrary irreducible moments of 
the distribution Pt„ (g) upon phonon excitation, as well 
as giving a picture of the underlying interplay between 
quasi-elastic and inelastic processes. 
The starting point of our derivation is the following 
model Hamiltonian 



H = ^H^, + ed^dld„+uj(,a'^a + V + Ve-ph (1) 
V" = ^ t^di^lkada + H.c ; Ve-ph = A(a + a^) ^ 

where a single molecular level of energy is coupled to 
the left (right) electrode by a hopping element iLd(ii?d), 
and interacts with a local phonon mode of energy 
with e-ph coupling constant A. The indexes (^, A:, a) la- 
bel the state of the ij, = L,R uncoupled electrode, char- 
acterized by wave vector k and spin a. We further de- 
fine the cumulant generating function (CGF) S{x) = 
~^n=i '''nl ('?")c ^ functional generating the ir- 
reducible moments of the distribution ('?")^- The con- 
nection of this definition to the former Hamiltonian is 
given by [1] 
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where V^^(t) denotes V with the substitution in the left 
hopping element tLd by tLrfe"'-^*-*-*/^, and Tc means time 
ordering on the Keldysh contour going forward from time 
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to time to and backward from time to to time 0. The 
counting field equals to ±x on the forward (back- 
ward) branch of the Keldysh contour and accounts for a 
virtual measurement of the charge being transmitted ■ 
As shown in Ref.[l3|, it is convenient to work with the 
generahzed current I{x) = ^i^'^^ix)^ that can be ex- 
pressed in terms of the Keldysh Green functions of the in- 
teracting molecular level G'^^{t,t') = —i{Tcd{ta)d'^ [t'^)) 

and of the uncoupled lead g^L (with indexes a, (3 = ±) 



27r 



I{X) = i^WTL j dj^[e}^G-+{u)gtL{oj) 

-'^G+-{u^)glt{u)] (3) 



In the former expression. 5 = 2 stands for spin degen- 
eracy and Tl — t\dl^ ^he coupling to the left con- 
tact expressed in units of the inverse of density of states 
W = l/irpL (supposed to be constant). 




At second order in perturbation due to e-ph interac- 
tion, the dot Green function can be written Gdd ~ 



(0) 
dd 



^m'^m^'^^m ^ ^^'^ the problem of finding the 
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CGF is thus equivalent to the one of computing e-ph 
self-energy in Keldysh space ff^d^^ the presence of 
the counting field. We retain two diagrams for the for- 
mer. The first one is a Hartree like diagram (FiglT](a)) 
which is diagonal in Keldysh space, frequency indepen- 
dent, of order A^/wqj and does not exhibit any jump at 
the inelastic threshold {V ~ ujo)- The second diagram 
is the exchange diagram (Fig[TJb)) that is responsible 
for the behaviour of transport properties at the inelastic 
threshold. The corresponding bubble expansion of S{x) 
is shown on FiglT](c) for the unperturbed CGF, FiglTfd) 
for the Hartree term and FiglTJc) for the exchange term. 
Taking into account the Keldysh indexes, one obtains a 
natural decomposition of I[x) as /o(x) + Iin{x) + Iei{x)i 
where Imix) is an inelastic contribution which arises 
from the non diagonal elements of the e-ph self-energy 
i^dd^ and E7^), and leiix) from the diagonal ones (Sj^ 



(b) 



and ^dtT) 1^1' This decomposition is equivalent to the 
one in Ref. 16 1. 

The unperturbed current /o(x) corresponds to resonant 
tunneling across the molecular junction in absence of e- 
ph interaction, and is given by 
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FIG. 1: Upper figure : Second order electron-phonon self en- 
ergy in Keldysh space (q,/3 — ±). (a) Hartree term, (b) 
Exchange term. Lower figure : corresponding bubble expan- 
sion of the CGF S{x)- (c) Unperturbed term, (d) Hartree 
term, (e) Exchange term. 
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{e'''.fLil-fR)-e-^^Ml-fL)} (4) 



+ {e^^ - 1)/l(1 - /r) + ie-'^ - l)/fl(l - U) 



s f dto 
2^ J 
1 
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where ,fL{B:) is the Fermi distribution of the left (right) 
lead, T{uj) = 4rLr_R/[r2 + Ed)^] the zero bias trans- 
mission coefficient and F = F^ +rfl the total coupling to 
the leads. The corresponding CGF coincides with the one 
derived by Levitov et al. [3|]. Effects of e-ph interaction 
are included in the two remaining terms. The inelastic 
contribution /i„ (x) can be written as 
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(e^'^F^/i + Fk/h)S-+ - [e-^^TUl - U) + Tn{l - /fl,))S+-] } 



(5) 



This corresponds to tunneling processes with absorption 
(emission) of a phonon. During such a process, the final 
energy of the scattered electrons increases (decreases) by 
an amount loo and the mean number of phonons decreases 
(increases) by one unit. The last term lei (x) accounts for 
elastic processes during which the energy of the scattered 



electrons is conserved, and is given by 



lelix) 



duj dA^ 
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UJ - UJd) 



y++ _ y — 

■^dd ^dd 



Fi(2/i-l) + Ffl(2/«-l) 



-'dd 



-'dd 



}(6) 
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The term involving T,'^^ — corresponds to a renor- 
malization of the transmission factor that gives logarith- 
mic corrections, whereas the term involving E^J*" + E^^" 
corresponds to quasi-elastic tunneling with emission- 
reabsorption of a phonon (hence the phonon population 
is unchanged) together with a virtual leaking of the prop- 
agating electrons into the leads. 

The former compact formulas can be used to implement 
a numerical computation of the FCS Of particu- 

lar interest is the behavior of I{x) at inelastic threshold, 
which is encoded in the jump of the generalized conduc- 
tance AG(x) = ^-f(x)L+ - At zero tem- 
perature (assuming phonon population /b =0), we find 
an analytical formula for AG{x) 



AG(x) 
AG,;„(x) 

AGe,(x) 



AG™(x) + AGei(x) 



(7) 
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2^ -e-P-r(g^)t 4r^r^ 
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Lr — L-i 
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e—ph r 

[(a-l)Li + (Q + l)(L3-i2)]} 



where we have introduced the dimcnsionlcss e-ph cou- 
pling Ag_p,j = A^/r^ and the parameter a = ^l/^r 
measuring the asymmetry in the coupling to the leads. 
The following notation is adopted A^.± = A^(±^), 
T± ^ T(±^), Li = 1/A^;+Ax;_{l/A^;+ + l/A^,_} 



and L 



2(3) 



This formula is the main result of the paper. We em- 
phasize that it encodes the full energy dependence of 
the model, i.e. it is not restricted to wide band ap- 
proximation, and allows to explore the behavior of the 
derivative with respect to voltage of the order n cumu- 
lant Fr, 



toT{o)\'i_^^ dv 
phonon energy given by 



AF„ 



(g")^, which exhibits a jump at 



AG(x) 



-ph 



9x""^ TiO) 
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We show in Figl2] the phase diagrams derived for the 
first three F„ factors, when exploring parameter space 
{a,T(0)} by modulating a = ^l/^r and shifting the 
molecular level e^. The first Fi factor exhibits a jump 
at phonon energy AFi (jump in the conductance), easily 
expressed as 



AFl;, 



T+T_ 



AFi 



l:el 



2Tr T(0)T(erf) 



27rT(0)T(ed) 2a 



} (9) 
T+} (10) 



The sign of Ai^i (represented on upper panel of Fig[21) 
has been studied in Ref.P,! [ill . [T^ . We find two regions 
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FIG. 2: (Color online) Sign of the jump in the derivative with 
respect to voltage of the order n cumulant at phonon energy 
AFn. Parameter space is {a,T(0)}, and phonon energy is 
Luo — lO^^r. From up to down : AFi (jump in the conduc- 
tance), AF2 (jump in the noise) and AF3 (jump in the third 
moment). The blue (red) color encodes a negative (positive) 
jump. 



of the parameter space corresponding to a negative jump 
(in blue) and a positive one (in red). If ^ T (which 
is the case for Figl2] where ujq = lO^^F), the contribu- 
tion due to inelastic processes AFi-m is positive when 
T{0) < 1/2, due to the opening of an inelastic channel 
and negative when T{0) > 1/2, due to enhanced inelastic 
backscattering (the correction to this strong coupling be- 
havior is of second order iiiujq/T). On the other hand, the 
contribution of quasi-elastic processes to the jump AFi-ei 
(of first order in wg/F) is always negative, due to elastic 
backscattering. Interestingly, AFi-ei is proportional to 
a — 1 and exactly zero when the contact is symmetric 
(a = 1). We emphasize that quasi-elastic processes con- 
tribute to the jump, because of virtual propagation into 
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the electrodes during the emission-reabsorption process, 
that is PauU blocked when V < luq (no available final 
scattering states). In the limit a ^ 0, both quasi-elastic 
and inelastic processes are of the same order of magni- 
tude and fully compete. The case of the second factor F2 
(middle panel of Figl2) corresponds to the jump in the 
noise at phonon energy AF2 



AF2:,. 



2{Tl + Tl 



T+T_)} 



■ {i-3(r+ 



27rT{0)T{ed\ 



a- 1 
2a 



[aT-{l - 2T- - T+) - T+(l - 2T+ 



F 



(11) 
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We find three regions of parameter space. For ujq <C 
r, and in the limiting case of a symmetric junction 
(a = 1), the resulting total jump in the noise AF2 = 
AF2-in + AF2-ei IS positive when T(0) < 1/2 - \/2^/2 or 
7^(0) > \/2 + \/2^/2 and negative otherwise. This change 
of sign can be understood by the following qualitative ar- 
guments. In the regime where T(0) goes to or 1, shot 
noise goes to zero due to Pauli principle, and activated 
e-ph interactions open an inelastic channel (positive con- 
tribution to noise). In the intermediate regime where 
T(0) « i, shot noise is maximal and activated e-ph in- 
teractions result in a negative contribution to noise. The 
same type of diagram is shown for the jump in the third 
moment Ai^3 (jump in the skewncss) on lower panel of 
Fig|21 where the competition between quasi-elastic and 
inelastic processes results in the partition of parameter 
space in four regions. 

The behavior of an arbitrary AFn and the whole FCS 
can be determined in the limit T ^ 1, where we obtain 
the following analytic approximation of the CGF [l^ 



S{x) ~ ~i%X~(li{e 



(13) 



where = toV/2T: and = to/2^\l_^f^/T{ed){V - 
llIq)9{V — uq). Below the inelastic threshold {V < ujq), 
the distribution Ptg (q) is Gaussian (a delta peak at zero 
temperature) due to perfect transfer of mean charge ^q, 
whereas above that threshold {V > ujq), Pto{q) is dis- 
torted to a Poisson distribution for holes due to the ac- 
tivation of sponteneous phonon emission (rare event for 
weak e-ph coupling). 

In conclusion, we have derived a compact formula for 
computing the FCS of a molecular level weakly inter- 
acting with a local phonon mode. The competition be- 
tween quasi-elastic and inelastic processes results in the 
partition of the parameter space {a,T(0)} into n + 1 
domains characterized by a given sign in the jump of 



the generalized cumulant of order n at phonon energy 
AFn- In the limit of perfect transmission, Ptoil) evolves 
from a Gaussian distribution for electrons to Poissonian 
distribution for holes, under activation of e-ph interac- 
tion. Of immediate experimental interest is the change 
of sign in the jump of noise. For temperatures in the 
range T = 4 — 17 K, and typical energy of the phonon 
mode ujo ~ 50 meV, the ratio T/wq « 10~^ - 4.10"^ is 
very small and the jump is not smeared by thermal ef- 
fects. The amplitude of the jump being of order a few 
percent, we expect that the change of sign could be tested 
experimentally along the lines of Ref . [ij] . 
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